Lieb-Liniger model of a dissipation-induced Tonks- Girardeau gas 
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We show that strong inelastic interactions between bosons in one dimension create a Tonks- 
Girardeau gas, much as in the case of elastic interactions. We derive a Markovian master equation 
that describes the loss caused by the inelastic collisions. This yields a loss rate equation and 
a dissipative Lieb-Liniger model for short times. We obtain an analytic expression for the pair 
correlation function in the limit of strong dissipation. Numerical calculations show how a diverging 
dissipation strength leads to a vanishing of the actual loss rate and renders an additional elastic 
part of the interaction irrelevant. 
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I. INTRODUCTION 

A Tonks-Girardeau gas is a one-dimensional (ID) 
system of identical bosons with strong interparticle in- 
teractions. The wave functions of the bosons have a 
surprisingly simple one-to-one mapping to a system of 
non-interacting fermions In recent years, a Tonks- 
Girardeau gas was realized in experiments with ultra- 
cold gases [3, 0] • A characteristic property of the Tonks- 
Girardeau gas is that the probability to find two parti- 
cles at the same position is very small Q . Theoretical 
and experimental studies of the Tonks-Girardeau gas so 
far dealt only with the case of elastic interactions. In 
a recent experiment, we studied a ID gas in which the 
bosons strongly interact inelastically, leading to loss of 
particles This system also shows a strong suppres- 
sion of the probability to find two particles at the same 
position and, indeed, demonstrates an alternative way of 
realizing a Tonks-Girardeau gas. The broader impact of 
this experiment lies in the demonstration that dissipation 
can be used to drive a system into the strongly-correlated 
regime. 

Here we study the theory of the dissipative Tonks- 
Girardeau gas. We show that in the limit of infinite 
interaction strength Girardeau's solution [2] is reached, 
as in the case of strong elastic interactions. We derive 
a Markovian master equation that describes the dissi- 
pation in the system and yields a rate equation for the 
particle loss. When the study is restricted to loss rates 
at short times, the master equation can be reduced to 
a Schrodinger equation with a non-Hermitian version of 
the Lieb-Liniger Hamiltonian Q. This approach yields 
an analytic expression for the loss rate in the limit of 
strong dissipation. Furthermore, we present numerical 
solutions to the dissipative Lieb-Liniger model at finite 
interaction strength that illustrate how a possible elastic 
part of the interactions becomes irrelevant as the dis- 
sipation strength diverges. In addition these numerical 
calculations reveal the dependence of the pair correla- 
tion function on dissipation strength, showing that an 
increase in the dissipation strength leads to a decrease of 
the actual loss rate due to the build-up of strong corre- 



lations. Finally, we present an estimate for the loss rate 
at longer times. 



II. INELASTIC ULTRACOLD COLLISIONS 

In this section, we discuss a master equation approach 
to describe loss processes caused by inelastic two-body 
collisions in an ultracold gas. First, we introduce the 
master equation, then use it to derive a rate equation for 
the particle loss, and finally show that in the short-time 
limit, a non-Hermitian Lieb-Liniger model is obtained. 



A. Master Equation 

We consider a three-dimensional (3D), dilute, ultracold 
gas of identical bosons, which all occupy the same inter- 
nal state. Two-body collisions in such a gas can be de- 
scribed in terms of the 3D s-wave scattering length a. For 
elastic collisions, it is customary to replace the interpar- 
ticle scattering potential a delta potential ^sd'^'^'^) (x) 
with interaction strength 



53D 



(1) 



where m is the mass of a particle. For the many-body 
system in second quantization, this yields a Hamiltonian 
(see appendix lAl or. e.g., Ref. [10]) 

Ho ^ J d3x*^(x)ff^(x)*(x) 

^Re(|3D) J rf3^^t2(^)q,2(^)^ (2) 

where i?s(x) = — ?i^V^/2m is the single-particle Hamil- 
tonian in first quantization and ^'(x) is the field operator 
that annihilates a boson at position x. In the case of 
purely elastic interactions g^u is real anyway, so taking 
the real part of it in Eq. ([2]) has no effect. 

We now generalize this treatment to include inelastic 
collisions. First, we consider inelastic collisions on the 
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level of a two-body problem. This yields an imaginary 
part of the scattering length, with Im(a) < [ll|, 
m, [l^. This is because the elastically scattered wave 
is proportional to the diagonal element of the S'-matrix 
g-2iafc„, ^ where krei is the wave vector in the relative 
motion. Hence, a negative imaginary part of a reduces 
the elastically scattered flux. This missing flux appears 
in the inelastically scattered channels. 

Next, we consider the effect of inelastic collisions on the 
level of the many-body system. The interaction strength 
f/SD defined by Eq. ([T]) becomes complex- valued. Equa- 
tion ^ contains only the real part of g^u and thus rep- 
resents only the elastic part of the interactions (see ap- 
pendix [A|. For a treatment of the inelastic collisions, 
it is crucial what happens to the collision partners after 
the collision. The difference in internal energy that is 
released in the change of the internal state appears in 
the kinetic energy of the relative motion of the particles 
after the collision. In typical experiments with trapped 
ultracold gases, the energy released here is so large that 
all particles involved in the collision quickly escape from 
the trap. This loss of particles is an irreversible process. 
In appendix \^ we derive a Markovian quantum master 
equation 



= -t[Ho,p]+Vip) 



(3) 



for the time evolution of the density matrix p of this 
many-body system, with a dissipator 



-Dip) = 



Im(53D) 



(4) 



that has a Lindblad form [l^, [l^ . 

The dissipator describes loss of pairs of particles due to 
inelastic collisions. More specifically, the last two terms 
in V deplete terms in the density matrix that represent 
states with more than one particle. The first term in V 
makes this lost population reappear in states with two 
fewer particles. This is further illustrated in Sec. IIIDI 

B. Loss Rate in 3D 

We now use this master equation approach to derive a 
rate equation for the loss of particle number. The particle 
density operator (x) = ^''l'(x)5'(x) has an expectation 
value that has a time dependence with a contribution 
from Hq , which conserves the total particle number, and 
a contribution from ^{p), which describes loss 



<i(n3D(x)) 



dt 



loss 



= -Tr[n3D(x)I?(p)] 

= -A'3D(*t2(x)*2(x)), 



where Tr denotes the trace and 
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-Im(53D) 
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(5) 
(6) 

(7) 



is the rate coefficient for two-body loss. The loss-rate 
equation ([6]) can be rewritten as 



d(n3D (x)) 



dt 



loss 



with the pair correlation function 



(^3d(x))2 



(8) 



(9) 



If N particles occupy the same single-particle wave func- 
tion, such as in a Bose-Einstein condensate (BEC), one 
obtains g^^^ — (N — 1)/^- For ^ oo such a sys- 



tem is uncorrelated, i.e., g 



(2) 



1. For comparison, a 



thermal gas of bosons above the critical temperature for 
Bose-Einstein condensation shows thermal bunching with 
g^^^ = 2. g^"^^ quantifies how far the loss rate deviates 
from that of an uncorrelated system. Similar considera- 
tions apply to inelastic three-body collisions HI, [3] ■ 



C. Loss Rate in ID 

The above master-equation approach, has a straight- 
forward generalization to a ID system, yielding 



d{hiD{x)) 



dt 



-K,u{niD{x))^g^^Hx) (10) 



loss 



Kid = -■^lm(giD), 



(11) 



where niD is the ID particle density and giD is the 
strength of the ID delta potential, which can be related 
to the 3D scattering length a as described now. 

Experimental realizations of ID systems typically use 
a strong harmonic confinement in the two transverse di- 
mensions, with trap angular frequency ujj_. If one approx- 
imates the 3D particle density along the two transverse 
dimensions by the ground state of the harmonic oscilla- 
tor with oscillator length aj_ = y^h/muj^, then spatial 
integration of Eq. ([8]) over the two transverse dimensions 
yields the estimate 



Km 



KsD 



(12) 



A more rigorous approach to this scatteriimproblem in 
a confined geometry was presented in Ref . [20] , yielding 



1 



\/2 aj 



C 



(13) 



where C denotes the Riemann zeta function with C(i) w 
— 1.46. In the limit \a\ ^ this reproduces Eq. il2p . 
But outside this regime, the scattering process is altered 
by the transverse confinement and Eq. (|12p is no longer 
valid. The derivation in Ref. considered only real- 
valued a, but following this derivation, one can show that 
Eq. ()13|1 remains valid for complex- valued a. The appli- 
cability of Eq. (|13p for the parameters of our experiment 
in Ref. [7] is discussed in appendix [Bl 
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D. Loss Rate at Short Times 

For studies of the behavior of the system at short times, 
it is useful to rewrite the master equation ([3]) as 



. dp 
I— 
dt 



iH.ffp + ipHlg + J{p) (14) 



with an effective Hamiltonian, that is not Hcrmitian and 
turns out to be the analytic continuation of Hq 



Im(53D) 



J d3x*t2(x)^r2(x) (15) 



and with 



J{p) = -Im(53D) y"d3a;^'2(x)p*t2(x). (16) 

Note that in a quantum Monte-Carlo approach to open 
quantum systems, J would be represented by the quan- 
tum jump part of the time evolution [2l[. 

The effect of the different terms in the master equation 
(|14p can be illustrated by considering a density matrix 
that initially represents a homogeneous BEC with exactly 
N particles. For short times Eq. (|14p yields 



p - (1 - Ct)\N) (N\ + Ct\N ~2){N~2\ + 0{et^), (17) 



where ^ = KsnnsoiN 
number is 



l)/2. The loss rate of particle 



(18) 



A detailed look at the calculation shows that Hcs 
causes the decay of the probability to find N particles, 
whereas ^7 causes the build-up of the probability to find 
N — 2 particles. If we are interested only in calculating 
dN/dt at short times, t <^ 1/^, we can use this obser- 
vation to drastically simplify the model. The key idea 
is that refilling population into states with N — 2 parti- 
cles is important only when the evolution at long times 
is considered, namely at times where states with N — 2 
particles also decay significantly. But at short times, all 
that counts is how fast population is lost from the initial 
state. 

Guided by this idea, we drastically simplify the model 
by dropping J' from the model. This removes the term 
^t\N-2){N-2\ fromEq. As a result, Tr(p) decays 

as a function of time, which is unphysical, but the rate 
at which this decay occurs 



dTijp) 
dt 



(19) 



is quite informative, because this is the rate at which 
population is lost from the initial state. We know that 
in reality this population reappears in states with N — 2 
particles. As each loss removes two particles, the true loss 
rate of particle number can be estimated from the decay 



rate of Tr(/3) in the simplified model. Indeed, comparison 
of Eqs. (UHl) and HH) yields 



dN 
'dt 



dTr(p) 



t=o 



dt 



(20) 



Note that here dN/dt refers to the true loss rate of par- 
ticle number in the full model, whereas dTv{p)/dt is a 
quantity calculated in the simplified model. The factor 
of 2 in this equation reflects the fact that two particles 
are lost per inelastic collision. 

With p initially prepared in a pure state p = 
and with J neglected, the master equation reduces to a 
Schrodinger equation with a non-Hermitian Hamiltonian 
Hcs- According to Eq. (|20|) the particle loss rate becomes 



dN 
~dt 



2 

ih 



en 



(21) 



For a right eigenvector of the effective Hamiltonian 

Hcs\'^) = E\'4)), we obtain 



dN 
~dt 



llm(i?). 



(22) 



The calculation of the initial loss rate is thus possible 
using a Schrodinger equation with iJcff, instead of the 
master equation. The imaginary part of an energy eigen- 
value E is always non-positive and represents the initial 
loss rate of the particle number. 

This greatly simplifies the calculation, because it suf- 
fices to study the eigenvectors and eigenvalues of an effec- 
tive Hamiltonian, instead of performing a time-resolved 
calculation of the density matrix. We will follow this 
approach throughout Sec. IIIIl 



III. DISSIPATIVE LIEB-LINIGER MODEL 

Motivated by the results of the previous section, we 
now study the eigenvectors and eigenvalues of the effec- 
tive Hamiltonian. This effective Hamiltonian is a non- 
Hermitian version of the Lieb-Liniger model We 
are mostly interested in the regime of strong dissipation. 
First, we show that in this regime, a Tonks-Girardeau gas 
is reached, second, we derive an analytical expression for 
the pair-correlation function in this regime and, third, 
we compare this with results of a numerical calculation, 
which show that a finite real part of the scattering length 
becomes irrelevant for diverging dissipation strength. 



A. Tonks-Girardeau Gas for Infinite Interaction 
Strength 

We consider a ID gas of A'^ identical bosons in a box 
of length L with periodic boundary conditions and with 



4 



the ID version of the effective Hamiltonian of Eq. p5l 
In first quantization this effective Hamiltonian reads 



2m ^ dxj 



giB ^ S{xt 



X,), (23) 



where Xi is the position of the ith boson. For the case of 
real- valued a, this model was first considered by Lieb and 
Liniger , who introduced the dimensionless interaction 
strength 



7 ^ 



(24) 



We explain now why this model in the limit I7I — > cxd 
yields exactly Girardeau's solution irrespective of 
whether the interactions are elastic or inelastic. 

The delta potential in Eq. (^5)) can be replaced by the 
following boundary condition at positions where the rel- 
ative coordinate Xij = xi — Xj vanishes Q 



dip 



dxj 



dip 



dxj 



mgiD 



Bosonic symmetry implies that on the left-hand side, the 
first term equals minus the second term. Expanding ip in 
a power series for 0, Eq. (^5)1 can be rewritten as 



i/'(xi, ...,xn) oc 



2r 

TO51D 



-fix,, I + 0(4), 



(26) 



where the proportionality contains the dependence on all 
remaining coordinates. 

The energy eigenstates of this model can be divided 
into two classes: gaseous states and bound states. As 
an example of a bound state, we consider a state where 
Re(7) < and where particles i and j are bound. Here 



ip oc exp 



/ ™giD | 



(27) 



for Xij 0, which obviously meets the condition (pS)) . 
An expansion in a basis of plane waves e*'^^'^ shows 
that this state contains imaginary momentum compo- 
nents that diverge for 7 — > — 00. Bound states exist only 
if Re(7) < 0. The case of real and negative 7 was studied, 
e.g., in Refs. [1, HI [H H [H, . 

For real 7, two characteristic properties can be used 
to distinguish between bound and gaseous states: First, 
the momenta of all particles are real valued for gaseous 
states, whereas at least one momentum is complex for 
bound states. Second, for 7 — s- ±00 all momenta converge 
to finite values for gaseous states, whereas at least one 
momentum diverges to ±ioo for bound states. 

We consider the case of complex 7, where all momenta 
are usually complex so that the first characteristic prop- 
erty becomes useless, but the second one remains useful, 
namely the convergence or divergence of the momenta for 
I7I ^ cxD can still be used to distinguish between gaseous 



and bound states. This distinction implies that ip and 
dip/dxij remain finite for I7I — > 00 for all gaseous states. 
Hence, the left-hand side of Eq. (P5)l remains finite for 
I7I ~^ 00, so that the right-hand side must remain finite, 
too. With giD diverging, the right-hand side implies that 



for 



I7l 



(28) 



for all gaseous states. The same result can be obtained 
from Eq. (|26p . because the term proportional to \xij\ 
must remain finite so that the proportionality factor must 
remain finite. With giD diverging, the constant term van- 
ishes. 

Yet another way to obtain Eq. is to consider a 
right eigenvector of the effective Hamiltonian i?cff IV'm) = 
Em\'>Pm)- If IgiB'tpi^ij = 0)| — > oo for at least one relative 
coordinate Xij, then Eq. (|23|) shows that Im(£') —00, 
so that the state decays with infinite speed. We decom- 
pose \Tp) — J2m'^m\''P'm) and sec that after any nonzero 
time only those eigenstates with finite gmipixij = 0) sur- 
vive. Again, this yields Eq. ((28)) . The other states are 
bound states. They vanish immediately and no longer 
contribute to the state or to the dynamics of the system. 

In the limit I7I —^ 00, the interaction for all gaseous 
states is fully described by the boundary condition pSjl . 
The crucial point is that this boundary condition is in- 
dependent of whether 7 is real or complex, and it is 
precisely this boundary condition that yields a Tonks- 
Girardeau gas. Hence, the wave functions of all gaseous 
states turn exactly into Girardeau's solutions 2] in the 
limit I7I — !■ 00, irrespective of whether 7 is real or com- 
plex. Attraction, repulsion, and dissipation all produce 
the Tonks-Girardeau gas in the limit of infinite interac- 
tion strength. 

The above discussion shows that an initially uncorre- 
lated state subject to strong dissipation will experience 
a rapid initial decay. Some fraction of the population 
will survive and this remaining population will be in the 
Tonks-Girardeau subspace, where the loss is slow. The 
loss will prevent the system from subsequently moving far 
away from the Tonks-Girardeau subspace. This freezing 
of the population can be understood intuitively with an 
analogy from classical optics [27| or as a manifestation 
of the quantum Zeno effect [3, [2^|. The focus of our 
study is the behavior of the system once it has reached 
the Tonks-Girardeau subspace, not the rapid initial de- 
cay. The experiment in Ref. [7[ started from a strongly 
correlated state already, so that no rapid initial decay 
occurred. 



B. Lieb-Liniger Solution 

Here, we summarize the central results of Lieb and 
Liniger 8] insofar as they are required for the following 
discussion. Derivations for all results of this section can 
be found in Ref. 8]. Following these derivations, one can 
show that these results also apply to the case of complex 
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D . The Lieb-Liniger solution is based on a Bethe ansatz 



iP{xi,...,xn) = a{P) exp 




,(29) 



which holds only for < a;i < 2:2 < ... < xn < L. 
For other values of the coordinates, the solution is ob- 
tained from bosonic symmetry and periodic boundary 
conditions. The sum extends over all permutations P 
of the numbers 1, N. The parameters ki, fcjv need 
to be determined. Once they are known, the amplitudes 
a{P) can be calculated easily 0]. The Bethe ansatz yields 
a nonzero solution for ip only if all kj are mutually dif- 
ferent [1] . The energy eigenvalue of this solution is 



^2 N 

E = ^y^ki 



(30) 



Lieb and Liniger characterize the solutions by quan- 
tum numbers 711, tiat^i, which are positive integers. 
For a given set of these quantum numbers the follow- 
ing set of coupled implicit equations is to be solved for 
j = l,...,iV-l [8] 



{kj+\ — kj)L — 27rrtj 



where 9ij abbreviates 



N 



7sj + l 



'1,3 



-2 arctan 



7niD 



), (31) 



(32) 



with |Re(6'ij)| < tt. This set of — 1 implicit equations 
determines only the A'^ — 1 quantities fc^+i — kj. There is 
one remaining degree of freedom. It can be expressed in 
terms of the quantum number Q 



no 



T ^ 



(33) 



which must also be an integer. For a given set of quantum 
numbers ni, ...,11^-1, the value of no is constrained by 
the additional condition 01 



If one restricts no to the range 



TV N 

— < nn < , 

2 " - 2 ' 




(34) 



(35) 



then there is exactly one value no for each combination 
of the ni, nAr_i. 

Adding N to no adds a momentum of 271 /L to each 
particle in the system. This corresponds to a center-of- 
mass motion of the g cLS 3(S Oil whole, which is trivial and of 



little interest. The values of no in the range ((35)l . how- 
ever, correspond to different internal quantum numbers 
and thus need to be taken into account 

While the Bethe ansatz does yield solutions for com- 
plex 7, it remains unclear whether a complete set of so- 
lutions is obtained. Our numerical results, discussed fur- 
ther below, suggest that each solution for real (positive 
or negative) 7 continuously deforms into exactly one so- 
lution for complex 7, thus suggesting that completeness 
for real 7 and complex 7 are closely linked. But a solid 
proof of completeness for complex 7 is beyond the scope 
of the present paper. On the other hand, the complete- 
ness issue has only limited relevance for the experiment 
discussed in Ref. [3], because here only the time evolu- 
tion of specifically-prepared gaseous initial states is to be 
tracked. 



Fermionized States in the Limit of Strong 
Interactions 



We now consider the limit I7I 00. We expand the 
arctan in Eq. (|32p in a power series and keep only the 
lowest order. This requires that (ki — kj)/nYD^ vanishes 
for large |7|. A sufficient condition is that all kj remain 
finite, which defines the gaseous states. 

For any large but finite I7I, the gaseous states can be 
divided into two classes: first, the so-called fermionized 
states for which \ {ki — fcj)/niD7| ^ 1 so that we can ex- 
pand the arctan and, second, the other states which must 
have very large kj so that they are highly excited. In the 
limit I7I — > 00 all gaseous states become fermionized. 

Insertion of this expansion into Eq. (|3ip yields 

- fc, = ^ ^n, + 0(|7r'). (36) 
L 7 -|- z 

Iterating this equation, we obtain for j = 2, N 



2 

L ^ + 2 ^ 



Insertion into Eq. ([33| yields 

AT-l 



(37) 



27r / 7 



k^ = no - 

7 + 2^ 

Insertion into Eq. ([34|) yields that 

1 ( ^^\.r . \ N-l 



0{\l\-') (39) 



must be an integer, thus fixing no modulo N . As no must 
generally be an integer, the term 0(171"^) must vanish 
automatically. For given values of 7 and ni, njv~i, we 
thus obtain no and all the kj . 
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D. Ground State in the Limit of Strong 
Interactions 



These general results for fermionized states can be used 
to calculate the ground state properties. The ground 
state of the gas phase is characterized by [1] 



rii = 71.2 = ... = riM-i = 1. 
Eqs. (1271) - ^ yield no = and 



kj — 



27r 



7 



i 7 + 2 



J 



N+1 



(40) 



(41) 



According to Eq. (jSO]) the ground state energy is 



E = N- 



6m 



7 + 2 



This result was previously derived by Lieb and Liniger 
Q in an alternative way for the case of real 7 and in 
the limit N 00. Our above calculation generalizes this 
result to complex 7 and finite N. 

Expansion to lowest order in 1 /j, taking the limit TV — > 
00, and insertion into Eq. ([^ yield the particle loss rate 
at short times 



dN 
'dt 



-N i2. ini 

3m 



oi\i\- 



(43) 



We use Im(l/7) = — Im(7)/|7p and compare with Eqs. 
(fTO)l and (fTTj) . This yields one of the central results of 
this paper 
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(2) 



3R 



(44) 



For the special case of real 7, this result was previously 
obtained in Ref. J^. Note that g^^^ for I7I 00, 
irrespective of whether 7 is real or complex. This is a 
consequence of the boundary condition ((28)) . 

Eq. (|15)) reveals that the loss rate depends on the scat- 
tering properties only in terms of Im(l/7) cx Im(l/5iD). 
We can rewrite Eq. (fT2)) as 



1 _ maj /I ^ C(l/2) \ 
5iD 2?i^ \a V2aj_) 



(45) 



and find that Im(l/giD) oc Im(l/a), irrespective of 
whether the condition |a| a_L is met or not. A possible 
alteration of the scattering process in the regime where 
the condition \a\ <C a± is not met affects giu in such a 
way that the loss rate in the fermionized regime is un- 
changed. In the fermionized regime, the dependence of 
the loss rate on the scattering length is always given by 



dN 
'dt 



cx Im ( — 

a 



(46) 



Note that in a case where an additional lattice is applied 
along the one dimension, the loss rate also obeys Eq. (|^S)) . 
see appendix IB] and Refs. (tI. [28|. 




FIG. 1: Real and imaginary parts of ki for the A'^ = 6 ground 
state with Re(7) = 0. The same data are plotted in a different 
way as the top- left curve in Fig. [2][a) . 



E. Numerical Solutions 

Solutions of the Lieb-Liniger model outside the limit 
I7I 00 can be found numerically. To this end, we 
numerically solve the set of coupled Eqs. ([5T|) and ([5^ 
with the standard root-finding algorithm in Mathemat- 
ica. Computation time issues restrict this approach to 
reasonably small N. It is not necessary to choose a spe- 
cific value of L for the numerical calculations, because 
the relevant equations can be written in terms of the di- 
mensionless parameters N , 7, and kL. 

We illustrate these solutions for the ground state char- 
acterized by Eq. (1^0)) with = 6 particles. Figure [1] 
shows the real and imaginary parts of one of the parame- 
ters, namely fci, as a function of — Im(7). An alternative 
way of displaying the same data is to show the trajec- 
tory that fci follows in the complex k plane when Im(7) 
is scanned from to —00. Such plots are shown in Fig. 
[5] for all kj . Parts (a) , (b) , and (c) each correspond to 
one fixed value of Re(7). For Im(7) — s- —00 the solutions 
in all parts of the figure converge to Girardeau's solution 
(shown as filled circles). The important result is that 
Re(7) becomes irrelevant for |Re(7)| <C |Im(7)| and that 
Girardeau's solution is reached in any case. 

The values of kj obtained in such numerical calcula- 
tions can be used to extract g^"^^ in the same way as in 
the derivation of Eq. (|^ . Results for the ground state 
are shown in Fig.[3Ka). g'^^^ depends on N, but converges 
in the thermodynamic limit. The solid line in Fig. [3Ka) 
shows the case N = 100. On the scale shown here it can- 
not be distinguished from the case A^ = 10 (not shown). 
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plane for the A'' = 6 ground state. Each curve shows a trajec- 
tory traced by kj when scanning the parameter Im(7) from 
to —oo. In all cases, the curves begin for Im(7) = some- 
where not too far from the coordinate origin and converge for 
Im(7) —oo to Girardeau's solution (shown as filled circles), 
(a) Re(7) =0. For 7 = all kj = 0. (b) Re(7) = 0.25. For 
Im(7) — 0, the kj are approximately equally spaced on the 



real axis, (c) Re(7) 



-0.25. For Im(7) = 0, the kj are 



approximately equally spaced on the imaginary axis, indicat- 
ing a bound state with an energy that is real and negative. 
However, as |Im(7)| grows and becomes larger than |Re(7)|, 
the finite Re(7) becomes irrelevant and Girardeau's solution 
is reached for I7I — > 00. 



This suggests that Fig. shows a reasonable approxi- 
mation to the thermodynamic hmit. A thorough analysis 
of finite-size effects is beyond the scope of this work. 

The loss rate Eq. (|10|) can be rewritten as 



dn 



ID 



2tinl 



ID 



dt 



Im(7)5(2). 



(47) 



The dependence on 7 is given by the dimensionless ex- 
pression Im(7)g(^). This quantity is displayed in Fig. 
[3Kb), showing that for — Im(7) 1 the loss suppresses 
itself. The experiment in Ref. Q was performed for 
I7I = 11 which is in this regime. 




FIG. 3: Pair correlation function and loss rate of the ground 
state as a function of the dissipation strength for Re(7) = 
and A'' — 100. (a) The pair correlation function g'^' (solid 
line) equals {N — 1)/N « 1 at 7 = and vanishes for 
7 — + ~ioo. For I7I ^ 1 it is well approximated by Eq. (I44p 
(dotted line), (b) The dimensionless quantity — Im(7)t;'-^-' is 
proportional to the particle loss rate, see Eq. (|47|l . For small 
I7I the loss rate increases linearly with Im(7) because g'^' is 
approximately constant. But for large I7I, the suppression of 



71 



dominates and the loss rate vanishes like 1/|7|. 



This is remarkable: the loss suppresses itself. 



IV. LOSS AT LONGER TIMES 

So far we studied only the loss rate at short times. A 
detailed study of loss rates at longer times would require 
a numerical solution of the master equation including 
the term J{p)- This is beyond the scope of the present 
paper, but for the case with a lattice applied along the 
one dimension, results of such a numerical calculation are 
reported in Ref. ^2^. 

Instead of solving the master equation, we now present 
a fairly simple estimate for the temporal evolution of the 
loss rate. During the loss, a spread in particle number 
builds up, but this spread is expected to reach a Poisson 
distribution, which is narrow enough that it suffices to 
consider only the mean value of the particle number. As 
the loss proceeds, the particle number changes and so 
does the ground-state energy Eg in Eq. (|42| . In the Lieb- 
Liniger model with periodic boundary conditions. Eg cx 
. In a more realistic model with a harmonic trap along 
the one dimension, the scaling is Eg oc N"^, which we 
use in the following. A simple estimate can be obtained 
when assuming that the average kinetic (i.e., real) part 
of the energy per particle is preserved by the loss, so 
that Jie{E)/N remains constant. We assume that E = 
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Eg at t = 0. As N decreases, Re{E)/Re{Eg) will then 
increase according to Re\E{t)]/Re[Eg{t)] = N{0)/N{t). 
The experiment in Ref. monitored a decrease of N 
by a factor of 2. Our simple estimate predicts an 
increase of Re[E{t)]/Re[Eg{t)] by the same factor of 2. 
We conclude that while the loss occurs, the energy of the 
system evolves somewhat away from the energy of the 
ground state, but not terribly far. 

This simple estimate is based on the ad hoc assumption 
that Re{E)/N is preserved by the loss. We now study 
how realistic this assumption is. As the loss takes the 
system into the regime Rc[_E(t)] > Re[Eg{t)], we need to 
study the excited states of the gas phase. We restrict our 
considerations to the fermionized regime. In a first step, 
we ignore the condition f34|) and set no = 0. Eqs. ([SO]) . 
dSZl), and dSHl) yield 



E\no=a — 



2m\L + 2 



f{ni,...,nN-i), (48) 



where 



N 



N-1 



j=l \s=l s=l / 

As the ground state truly has uq = we obtain 
/(ni, ...,nN-i) 



E\ 



no=0 



En 



(49) 



(50) 



where Eg denotes the ground-state energy. Obviously / 
is always real. Hence, 



Im(£;|„„^o) ^ ImjEg) 
Rc(S|„„=o) RciEgY 



(51) 



We will show now that the correction for no 7^ is neg- 
ligible. To this end, we note that according to Eqs. ([57|) 
and ([55]) . no shifts all kj by 2TrnQ/NL. As a result, the 
energy becomes 



E — E\n„=o + 



1%^ / 27rno \ ' 
2mN \ L ) 



(52) 



Eq. ([35]) yields |no|/i < niD/2. In the thermodynamic 
limit {N ^ 00 at fixed niD), the correction is of order 
E-E\na=o <nl^O{l/N). On the other hand, i;|„o=o > 
Eg cx NuIy). Hence, 



E 



E 



no=0 



1 + 0{N- 



(53) 



Therefore the correction due to no is negligible and Eq. 
(l5T]l becomes 



Im(£') _ Im(£;g) 
Re{E) ~ Re{Eg)' 



(54) 



This shows that the ratio of the real to the imaginary part 
of E is identical for all fermionized states, independent 



of ni, ...,n7v_i. Combined with Eq. (|22|) we obtain that 
the loss rate of a state is proportional to the kinetic {i.e., 
real) part of the energy of a state. 

This has two important consequences for the under- 
standing of the loss at longer times. First, it shows that 
an increase of Re[E{t)]/Re[Eg{t)] by a given factor leads 
to an increase of the loss rate compared to that of the 
ground state by the same factor. Second, it shows that 
the loss preferentially depletes states with large energy. 
Hence, the above estimate assuming that the loss leaves 
the average kinetic energy per particle unchanged is too 
pessimistic. In reality, the ratio Re[E{t)]/Re[Eg{t)] in- 
creases even more slowly than predicted by the simple 
estimate above. The kinetic energy and loss rate of the 
system will thus remain fairly close to those of the ground 
state if the loss is allowed to evolve only for a time where 
the particle number decays by a factor of ~ 2 as in Ref. 



V. CONCLUSION 

In conclusion, we showed that strong inelastic interac- 
tions produce a Tonks-Girardeau gas, much as strong 
elastic interactions. We derived an analytic expres- 
sion for the pair correlation function of the fermionized 
ground state. Numerical results illustrate the behavior 
of the ground state for small and medium dissipation 
strength. These calculations also show that the loss sup- 
presses itself. Finally, we presented a simple estimate for 
the evolution of the loss rate for longer times. 



APPENDIX A: DERIVATION OF THE MASTER 
EQUATION 

1. General Considerations 

In this appendix we derive the master equation ([3]) 
with the dissipator ([4]) for the description of loss of par- 
ticles from a many-body system due to ultracold inelas- 
tic two-body collisions. The system of interest S consists 
of trapped motional states in the initial internal state. 
Inelastic collisions couple this system S to other states, 
which we call the reservoir R. The reservoir states can 
represent all possible internal states. The motional parts 
of the reservoir states are untrapped because we assumed 
that all particles involved in an inelastic collision quickly 
escape from the trap. Hence the energy spectrum of these 
untrapped states is a continuum. 

The total system consisting of S and R is represented 
by a density matrix x. As we are not interested in the 
reservoir, we can trace out the reservoir degrees of free- 
dom and obtain a density matrix p that describes only 
the system of interest S. We assume that the reservoir is 
in the vacuum state initially and we use a Born-Markov 
approximation for the reservoir modes, i.e., we as- 
sume that the reservoir has no memory of the past. This 
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is justified because the reservoir consists of a continuum 
of modes that quickly dephase. As a result, the reser- 
voir returns to the vacuum state on a time scale, which 
is much faster than all time scales in the system of in- 
terest S. Using very general arguments one can show 
[isl that such a Born-Markov approximation yields a 
Markovian quantum master equation 



^1 ^ -[-^^ 



(Al) 



for the evolution of the density matrix p describing only 
the system of interest S. This master equation contains 
a Hamiltonian Hi and a so-called dissipator P, which 
models the irreversible part of the time evolution and 
can always be written in a Lindblad form [l^ [l^ 



AlAup- 



pAlAk) , (A2) 



where the 7^ > are relaxation rate coefficients and the 
Ak are called Lindblad operators. 

The rest of this appendix is devoted to deriving ex- 
pressions for Hi and 2? for our particular system. Given 
the form of the elastic collision term in Eq. the fi- 
nal result for the Lindblad operator A^ = ^'^(x) in Eq. 
(HI is not terribly surprising. But the result for the rate 
coefficient jk = — Ini(53D)/^ is certainly less obvious, in 
particular as far as possible factors of 2 are concerned. 
We now present a rigorous derivation for this dissipator, 
including the rate coefficient. The calculation proceeds 
in three steps. First, we determine a contact potential 
that can be used to describe inelastic scattering, second, 
we derive the corresponding Hamiltonian in second quan- 
tization and, third, we eliminate the reservoir degrees of 
freedom. 



the relative motion hats and lia'fS' , respectively. The wave 
function (|A3D can be regarded as a generalization of Eq. 
(5) in Ref. 29] to inelastic s-wave scattering. The fact 
that the outgoing particle flux is proportional to ka/pi 
makes it necessary to include the factor ^/kap/ka'p' in 
Eq. (|A3[) , in order to maintain the fact that the 5'-matrix 
is unitary if the number of particles is conserved in the 
collision. 

We now replace the interparticle interaction potential 
by a contact potential. The wave function (|A3P that is 
usually valid for r — > cxd only, then becomes valid for all 
r. We choose a contact potential V^^^^ = gS^V^ with a part 

g that acts only on the internal state and a part (5icg that 
acts only on the spatial part of the wave function. For the 
latter, we use the regularized delta function (r|5rcg|^) — 
(5(3)(r)|:(rV') (see, e.g., Refs. [1, H^). In addition, we 
assume that g is Hermitian. 

We expand the wave function (|A3|) in a power series for 
r — !■ and spatially integrate the Schrodinger equation 
(r|(_ff — = over a sphere centered around r = 

with radius e — > 0+. We use V^i = ~Att5'^'^\v) and 

drAiTr'^5^^\Y) — 1. We neglect terms 0(e). This yields 
for aU aPa'(3' 

ga' fi' ,a" 13" ^/ ka" 13" {Sa" j3" ^af} + Sa"aSf3"l3) = 

a"l3" 



27rn^ 



IJ-yka 



--{Sa'f3',af3 — Sa'a^p'/s) 



(A4) 



with the reduced mass p = m/2 and the Kronecker sym- 
bol S. We do not attempt to solve this linear system for 
the variables ga'p',af3 = {<^' (3'\g\ctf3) in full generality. In- 
stead, we decide to fulfill Eq. ()A4|) for only one specific 
initial state \a(3) — We thus assume 



2. Contact Potential for Inelastic Scattering 

For the first step, we consider a scattering process with 
two distinguishable incoming particles and assume, for 
simplicity, that after the collision, there are two outgoing 
particles with the same mass as the incoming particles. 
We further assume that the incoming wave is an s wave 
and that the potential is spherically symmetric, so that 
the outgoing wave is also an s wave. The scattering state 
IV') at large relative distance of the particles is 




■Ja'l3',al3 — 



\a'f3'), (A3) 



where r is the relative coordinate, Sa'f3\ap is a matrix 
element of the scattering matrix, and the internal states 
of the particles before and after the collision are jo;/?) = 
\a) (8) 1/3) and |a'/3') with corresponding wave vectors for 



ga'P',ai3 = if (a',/3') and (a,/3) 7^ (i,i). (A5) 

This yields 



p{Sii^ii + 1) kapki', 
gii,a[3 — gaf3,ii 



gi; 



pkii \ Sii^ii -t- 1 p 



(A7) 
(A8) 



where we used that unitarity of the S'-matrix implies 

The scattering phase rj, defined by Su.a — e^*'', can be 
used to rewrite the above results in the compact form 



2nh Re(tan77) 
ga.ii = 7 (A9) 

p Kii 

E - ('-^)''^^. (AlO) 



^11 
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For ultracold collisions, the ratio tanrj/ku becomes con- 
stant and is called the scattering length 



a = — lim 



tan 77 



(All) 



The above treatment applies to scattering of distin- 
guishable particles, while what we are really interested in 
is scattering of identical bosons. Due to the linearity of 
the Schrodinger equation, scattering of identical bosons 
is closely related the scattering of distinguishable parti- 
cles from the same potential. Our following approach is 
based on second quantization, which has bosonic symme- 
try built into it automatically, so that no separate treat- 
ment of scattering of identical bosons is needed here. 



d3a:^-1^,(x)^'t (x)*„(x)^'^(x), (A14) 



where the field operator 5'3j(x) — J^k'^akiM^) creates 
a boson at position x in internal state \a). 

Note that the term in Hsc that describes elastic scat- 
tering in the initial internal state is proportional to gu^u 
from Eq. (jA9|) . Using Eq. (IAll|) . we obtain exactly the 
elastic collision term in Eq. ((2]). The resulting mean- field 
energy is small compared to the kinetic energies involved 
in the following calculation. Hence we can perform an ap- 
proximation of independent rates of variation (3^ . This 
means that we include the elastic collision term in the 
final Eq. ([2]), but we neglect its effect on V by setting 
gii,ii — in the following calculation. 



3. Second Quantization 

The interactions described by the above contact poten- 
tial can be translated to a Hamiltonian in second quan- 
tization for a many-body system. If a general two-body 
interaction potential y*^^) (x^, x^) is summed over all pos- 
sible pairs of particles, one obtains X];i<i/ ^^^H^M'^i^) ~ 
^ ^''^H^M' ^!^)' where /i, v enumerate the particles. 

In second quantization, this yields a Hamiltonian [3ll | 

Hv = IT. {j,k\V^'^\l,rn)blblbib^, (A12) 

jklm 

where the states \j) form a basis of single-particle states 
and 6] creates a boson in state 

We now apply this to our contact potential. As a ba- 
sis of single particle states, we choose |a) |k), where 
I a) is an internal state and |k) a momentum state with 
(x|k) = g«k x^-3/2 in box quantization with quantization 
volume L^. We denote the corresponding creation oper- 
ator as ajjj.. A short calculation yields the scattering 
Hamiltonian for our contact potential 

Tj _ \ ' \ ' \ ' 9a'l3'.afi t t 

X apM-v\ dVe-»(''+i) "-43)(r)e«(k+p) r^(A13) 



where 2?ik is the center-of-mass momentum of a pair of 
bosons, and Tip and Tiq are the relative momenta before 
and after the collision. We denote the corresponding ki- 
netic energies as — {2hk)'^ /4m, hujp = Ti'p^ /2n, and 
TiLOq = Ti'f^ Note that in Eq. (|A13p it makes no 
difference if we take 5^ol(r) or 5(3) (j.). 

Let us assume for a moment that the ga'[3',ai3 in- 
dependent of k, p, and q. Hgc can then be transformed 
into a simple expression in position representation 



4. Eliminating the Reservoir 

Having determined the explicit form of i/sc, we now 
proceed to eliminate the reservoir modes. To this end, 
we split the total Hamiltonian of our problem into three 
parts: Hs acting only on the system of interests S*, Hn 
acting only on the reservoir i?, and i?s_R coupling S and 
R. The single-particle Hamiltonian is thus decomposed 
into 

Hs+Hr = ^{— — + j ai^kOa.k, (A15) 

where the term a = i belongs to Hg, all others to Hji. 
Here TiWq is the internal energy in internal state |a). We 
choose the zero of internal energies such that oji = 0. 

For simplicity, we neglect collisions except if they make 
a transition between S and R. Hence Hsr = Hsc with 
the ga'i3',ap from Eqs. (|A5|) - (|A7|) and with the assump- 
tion ga^ii = 0. We further simplify the model by ex- 
cluding the possibility that a reservoir state contains any 
particles in the initial internal state This is a simple 
way to reflect the fact that all particles involved in an 
inelastic collision are lost from the sample. 

Our discussion throughout the rest of this section 
closely follows chapter 1 of Ref. [lH|. The total sys- 
tem consisting of S and R is described by a density ma- 
trix that obeys the von-Neumann equation ifix — 
[Hs + Hr + Hsr, x]- We transform this into an interac- 
tion picture 

X{t) = e^(^^+^«)*/''x(i)e"^^''^~+''«^*/^ (A16) 
HsR{t) - e^(«-+^«)*/''i7sfle-'(^^+^«)*/'' (A17) 

^ = ^[HsR{t),m)]- (A18) 
iTi 

Note that Hsr picks up an explicit time dependence. We 
formally integrate the last equation 

m = m + ^f^dt'[HsR{t'),x{t')] (A19) 
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and substitute this result into the commutator in Eq. 

(Eni) 

in 

_J_ [ dt'[HsR{t),[HsR{t'),x{t')]]. (A20) 
n Jo 

We assume that R is in the vacuum state at i = 0, uncor- 
related with the state of S. We assume that Hsr is weak 
and that the reservoir R is much larger than the system 
S, so that the state of R is hardly affected by the inter- 
action with 5*. The correlations between S and R that 
build up during the time evolution are then weak. In 
a Born approximation ^151] , we neglect these correlations 
and write x{t) — p{t)Roi where Rq — |0)(0| is the vacuum 
density matrix for R. We insert this into the right-hand 
side of the last equation. We then take the partial trace 
Trfl over the reservoir, assume Tvii[HsR{i)Ro\ = and 
obtain llal 



P = 'T2 I dt'TvR{[HsR{t),[HsR{t')rp{t')Ro\]}.{A2l) 
n Jo 



We write the scattering Hamiltonian (|A13[) with the 
ga'i3',ap from Eqs. (|A5|) ~ (|A7|) and with the assumption 

gii,ii = as 



Hsr = hJ2^lsi, + H.c. 



— 1i,k+pfli,k-p 



(A22) 
(A23) 



^k = EE^"/3«i,k+q4,k-q (A24) 
a/3 q 



(A25) 



We recall our assumption that a reservoir state cannot 
have any particles in state \i). Hence, Sk acts only on 
S (internal state whereas Fk acts only on R (all 
other internal states). In addition, our above assumption 
Tj:ji[HsR{t)Ro] = is satisfied. Moreover, ga'p',ap = 
gi3'a',i3a bccause swapping which internal state is labeled 
first, cannot change the physics. Hence Kap = Kpa- 

Insertion of Hsr into Eq. (|A21[) shows that Tr/? van- 
ishes for most of the terms. We obtain 

P^- f rft' V[sl,(t),5k(Op(i')]fik(i,t') + ^-c.,(A26) 

where we abbreviated Rk{t,t') = Tr;^{fk(i)f k(t')^o}- 
We use aa.},{t) = ,^e-*('^"+'"'*='/2")* and obtain 

Rk{t,t') = E E E 

Q/3 a' f3' qq' 

^ ^-i{LOa,+Ull3+Ulk+^q)t + i{LU^,+UIi^,+LUk+'^ql)t' 

X (0|aa,k+qa/3,k-qal.,k+q'4',k-q'|0>' (A27) 



For a 7^ /3 or q 7^ (or both) the last line equals 

(A28) 

Using Kaf} = KfSa wc scc that both terms produce iden- 
tical results after carrying out X^a'/s'q'- -f'or a = P and 
q = 0, we obtain only one term, but after carrying out 
Sa'/j'q'' the result is the same because (0|a^at^|0) ~ 2. 
Hence, 



(A29) 



a/3 



We now employ a continuum approximation for q. At 
a given center-of-mass momentum, the density of states 
in q-space is {L/2Tr)^. To see this, consider first a 
ID system. The possible values of the single-particle 
wave vectors ki or k2 are separated by 27r/L. Assume 
that a pair (fcijfe) corresponds to a certain center-of- 
mass momentum. The next possible pair of values that 
corresponds to the same center-of-mass momentum is 
{ki +L/2Tr, k2 — L/2tt), which changes q = {ki —k2)/2 by 
L/27r. Conversion into 3D and then into frequency space 
yields the density of states 



L\ A 2 P 

2^) ^"'^ ^' 



(A30) 



where fi/hq — dq/dtOq. Hence 
i?k(i,t') 



2 gi(Wa+C/J/J+Wfc)(t'-t) 



a/3 



X / du;qGiu;q)\K^0\^e'^^^''~'\ (A31) 



G{ujq)\Kaf3\'^ is typically a slowly varying function of ujq. 
For large t — t', the phase factor e*"'*^* oscillates 
rapidly and thus the integral is almost zero. This is the 
rapid dephasing of the reservoir stats that we mentioned 
earlier. The integral over t' in Eq. (jA26P is thus dom- 
inated by the values at t' « t. This justifies a Markov 
approximation [15], which consists of replacing p(t') in 
Eq. (|A26p by p{t). Since the integral over t' is domi- 
nated by the values at w i, we can extend the lower 
bound of the integral to t' = — oo [IB]- We substitute 
T — t~t' and Eq. (|A26|) becomes 

P - E [4k+p4k-p'a^k+p'aj,k-p'P(i)]Bkpp'(i) 

kpp' 

+H.C. (A32) 



with 
■Skpp' {t) 



drRuit, t')e'('^*=+'^p)*-»('^*.+'^p')t' (A33) 



Q/3 

xG{ujq)\Ka(3\'^e' 



dT 



duJn 



(A34) 
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We swap the two integrals, performing the time integra- 
tion first. We use llSll 



dre 



TrS(Lu) 



(A35) 



where P denotes the principal value for the following duj 
integration. Hence 



Id 
2 



' 2 



= 27r^ G(tJ,)|Ka 



/3| 



a/3 



a/3 



G'(^g)|«a;3pcg^g 



(A36) 
(A37) 

,(A38) 



where G{LOq)\Kap\^ in Eq. (|A37|) is to be taken at = 
Ljpi — ~ ^^,31 which reflects energy conservation. 

Note that our results for 7D and A are twice as large 
as Eqs. (1.68) and (1.70) in Ref_^. This difference is 
due to the two terms in Eq. (jA28[) . Ref. [lB| considers 
reservoir states with only one particle, so that only one 
term occurs in the expression analogous to Eq. (jA28p . 

We transform Eq. (jA32|l back from the interaction pic- 
ture to the Schrodinger picture. We use that 7^) and A 
are real and obtain 



^ 2 

kpp' 



k+p'^i:,k-p 



with 



kpp' 



,k+p'*i,k-p'^i,k+p'aj,k-p'- 



(A39) 
H.c. 

(A40) 



H/^ is very reminiscent of the energy shift in ordinary 
second-order perturbation theory, because H/^ describes 
processes in which a pair of particles in state \ii) makes 
a transition to a state \a(3) with relative momentum hq 
and then returns to state Energy is usually not con- 
served in such virtual transitions and the energy differ- 
ence between the initial and intermediate state appears 
in the denominator of Eq. (|A38[) . This denominator gives 
the typical weighting for such virtual transitions, known 
from perturbation theory. 

A calculation of the value of A from Eq. (|A38[) is dif- 
ficult, even if the incoming particles are ultracold. This 
is because, first, Eq. (jAlOp cannot be used because the 
energy denominator depends on a and f3 and, second, 
because energy is not conserved in the virtual transitions 
so that a calculation of the integral requires knowledge 
of the ^-matrix elements in Eq. (|A6[) over a large range 
of energies hujq, even if tkop' is fixed to be ultracold. The 
energy dependence of the ^-matrix elements is sensitive 
to many details of the true potential. A calculation of A 
is thus beyond the scope of the present work. 



We see that Ha in Eq. (|A40|) has exactly the same form 
as the elastic collision term {a, (3) = (a', (3') = (i, i) in Eq. 
(|A13|) . In the ultracold limit, thus simply contributes 
to the mean-field energy. Any attempt to measure Re(a) 
in the many-body system will always yield the sum of 
the two contributions, so that we can simply absorb A in 
an effective value of Re(a), that needs to be determined 
exp er iment ally. 

A calculation of from Eq. (jA37[) is much simpler 
because it requires knowledge of the S'-matrix elements 
only for outgoing momenta Ti(\ that satisfy energy con- 
servation. There is no additional denominator and we 
can use Eqs. jSni, and ((X30)) with q = k^/j and 

p' = kii. We obtain 



Id 



Airh Im(tan7y) 



(A41) 



For ultracold incoming particles, we use Eqs. (jAllj) and 
^ and obtain 



Id 



Ini(ff3D) 



(A42) 



This is independent of k, p, and q, so that Eq. (|A39|) 
can be transformed into a simple expression in position 
representation. This yields Eq. which is the central 
result of this appendix. 



5. Comparison with Existing Literature 

Finally, we discuss how our results relate to the existing 
literature on scattering theory. To this end, we note that 
the rate coefficient ifsD in Eq. ([7]) can be calculated in 
an alternative way, based on the observation that the 
scattering rate per particle in a gas is proportional to 
n^Yyahkre\/m, where a is the cross section and kj-ei is the 
wave vector of the relative motion of the two particles 
before the collision. For elastic and inelastic two-body 
collisions of identical bosons one obtains [13, uM 



CTci = 87r|a 

finel = 



— Im(a), 



(A43) 
(A44) 



respectively. Deriving the rate coefficient ii'sD from fiinci 
is a delicate issue because of several subtle factors of 2. 
An unambiguous way to handle those factors is based 
on a quantum Boltzmann equation (33j and confirms our 
resuh Eq. 



APPENDIX B: APPLICABILITY OF THE 
MODEL TO THE EXPERIMENT IN REF. 

This appendix discusses the applicability of the model 
Dresented in Sec. Ill Cl to our experiment reported in Ref. 
3. First, we note that the derivation of gxo in Ref. |23| 
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relies on the assumption that the kinetic energy in the rel- 
ative motion of two particles h?k'^^^/m before a colhsions 
is less than the splitting of the transverse harmonic os- 
cillator levels with even angular momentum 2hLo^. This 
is equivalent to fcrcia± < V^. 

In a fermionized ground state, noticeable momenta can 
occur even at zero temperature. We can use the kj to 
estimate the typical momenta in the ground state. We 
note that the kj are not actually momenta, see e.g. Ref. 
(3^ , but the typical width of the momentum distribution 
is set by kj. We are interested in the typical relative 
momentum in the fermionized ground state, which we 
can thus estimate from Eq. (|37p as fcroi ^ /cat — fci — 
{N — 1)2t: / L w 27rniD. The experiment in Ref. 0] was 
performed at an initial density of niD = 2/(830 nm) and 
with aj_ = 39 nm, so that fcj-eiaj. ^ 0.60 which is well 

inside the regime fcrcia± < \/2. 

Second, we mention another assumption in the deriva- 
tion of giD in Ref. 20], namely that the ID scattering 
process can be modeled by a ID delta potential only for 
low enough momenta. For fcreioj. ^ 0.60, Re(a) — 0, and 
K^D — 1.5 X 10^^" cm'^ I s 1351 we numerically evaluate the 
analytic results of Ref. [23| . We find that the transmis- 
sion coefficient [2^ for the scattering process is altered 
by less than 10% when making the transition from the 
3D regularized delta potential to the ID delta potential, 
thus showing that the model is, indeed, a good approxi- 
mation. 

Finally, we briefly discuss the part of our experiment 
in Ref. J\ in which a weak lattice is applied along the 
one dimension. If two particles occupy the same lattice 
site, the resulting on-site interaction is described by a 



matrix element U . We will discuss now, why the ex- 
pression U = gsn J (Px |w(x)|'' is valid, where u;(x) is a 
Wannier function. First, the lattice potential at a single 
lattice site is to a good approximation harmonic, with an 
aspect ratio of less than 3 for the parameters of our ex- 
periment. Scattering processes within one lattice site are 
therefore not in a ID regime and only for much larger 
aspect ratios would we obtain U = giD J dx \w{x)\'^, 
with a ID Wannier function. Second, if two particles 
occupy the same site, this 3D system at one site has 
|ap?^3£) <C 1. To see this we consider a lattice depth of 

Vo = mSr 

along all three dimensions, where Er is 
the recoil energy, yielding a harmonic oscillator length 
of aho = 39 nm at one site. At K^d = 1.5 x 10^^° 
cm'^/s [3^ and Re(a) = 0, we obtain for two particles 
at one site \a\^n3D,pciik = 2\a\^ / {^/iraho)^ = 0.03. The 
fact that |apn3£) <C 1 means that two particles at an 
isolated single site are in the weakly correlated regime 
and the single-particle Wannier function can be used to 
calculate U. 
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